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The Monte Carlo HYDJETH — h model, that combines parametrized hydrodynamics with jets, is 
employed to study formation of second V2 and fourth V4 components of the anisotropic flow in 
ultrarelativistic heavy-ion collisions at energies of the Relativistic Heavy Ion Collider (RHIC) and 
the Large Hadron Collider (LHC), -/s = 2QQA GeV and ^ = 2.76^ TeV, respectively. It is 
shown that the quenched jets contribute to the soft part of the V2{pt) and 114 (pt) spectra. The jets 
increase the ratio V4/V2 thus leading to deviations of the ratio from the value of 0.5 predicted by 
the ideal hydrodynamics. Together with the event-by-event fluctuations, the influence of jets can 
explain quantitatively the ratio V4/V2 at pr < 2GeV/c for both energies and qualitatively the rise 
of its high-pT tail at LHC. Jets are also responsible for violation of the number-of-constituent-quark 
(NCQ) scaling at LHC despite the fact that the scaling is fulfilled for the hydro- part of particle 
spectra. 

PACS numbers; 25.75.-q, 25.75.Ld, 24.10.Nz, 25.75.Bh 
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I. INTRODUCTION 

The transverse collective flow of particles is an impor- 
tant characteristic of ultrarelativistic heavy-ion collisions 
because the flow is able to carry information about the 
early stage of the reaction. Particularly, the collective 
flow is very sensitive to change of the equation of state 
(EOS), e.g., during the quark-hadron phase transition. 
The azimuthal distribution of particles can be cast [1, 2] 
in the form of Fourier series 
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Here pt, and y are the azimuthal angle, the trans- 
verse momentum, and the rapidity of a particle, re- 
spectively. The unity in the parentheses represents the 
isotropic radial flow, whereas the sum of harmonics refers 
to anisotropic flow. The first two harmonics of the 
anisotropic flow, dubbed directed flow vi and elliptic flow 
V2^ have been extensively studied both experimentally 
and theoretically in the last 15 years (see, e.g., [3] and 
references therein), while the systematic study of higher 
harmonics began quite recently [4-7]. 

In the present paper we investigate the ratio R = W4/f | 
in heavy-ion collisions at energies of the Relativistic 
Heavy Ion Colhder (RHIC) (^s = 200^ GeV) and the 
Large Hadron Collider (LHC) (Vs = 2.76yl TeV). Inter- 
est in the study was raised due to the obvious discrepancy 
between the theoretical estimates and the experimental 
measurements. On the one hand, the exact theoretical 
result for hydrodynamics provided V4/V2 = 0.5 for a ther- 
mal freeze-out distribution [8]. On the other hand, it 
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was found soon in RHIC experiments [9, 10] that the 
measured ratio R exceeded by factor 2 the theoretically 
predicted one. Both the STAR and the PHENIX Collab- 
orations have reported that the R is rather close to unity 
for all identifled particles in a broad ranges of central- 
ity, 10% < ajogeo < 70%, and transverse momentum, 
Pt > 0.5GeV/c. For the smaller pT the ratio seems to 
exceed the value of 1. Note also that the PHENIX data 
are about 10-15% below the STAR ones. 

In Ref. [11] it was argued that the experimentally mea- 
sured R can be larger than 0.5 even if the ratio V4/V2 
was exactly equal to 0.5 in each event. Such a distortion 
can be caused by event-by-event fluctuations. Namely, 
if the ratio V4/V2 is estimated not on an event-by-event 
basis but rather on averaging of both V2 and V4 over the 
whole statistics, the event-by-event fluctuations will sig- 
nificantly increase the extracted value of the ratio. Cal- 
culations of R at RHIC energies within both ideal and 
viscous hydrodynamics with different initial conditions 
[12] revealed that the ideal hydrodynamics provided bet- 
ter agreement with the data, although the STAR results 
remained underpredicted a bit. For LHC the hydrody- 
namic calculations have predicted similar behavior with 
slight increase at small transverse momenta [12]. 

The preliminary results obtained in Pb -I- Pb collisions 
at y/s — 2.76A TeV favor further increase of the V4/V2 
ratio [5, 6]. Moreover, this ratio is not a constant at 
Pt > 0.5GeV/c but increases with rising transverse mo- 
mentum. The first aim of the present paper is to study 
to what extent the hard processes, i.e., jets, can affect 
the ratio R predicted by the hydrodynamic calculations. 

The second aim of the paper is investigation of the 
fulfillment of the so-called number-of-constituent-quark 
(NCQ) scaling, observed initially for the partial elliptic 
of mesons and baryons at RHIC [13, 14]. Despite the gen- 
eral expectations, the measurements show that the NCQ 
scaling is broken at LHC energies [15]. Thus, it would 
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be interesting to elucidate the role of jets in the scaling 
violation. For these purposes we employ the HYDJET++ 
model [16], which couples the parametrized hydrodynam- 
ics to jets. The soft part of the HYDJET++ simulated 
event represents the thermalized hadronic state where 
particle multiplicities are determined under assumption 
of thermal equilibrium. Hadrons are produced on the hy- 
persurface, represented by a parametrization of relativis- 
tic hydrodynamics with given freeze-out conditions. At 
the freeze-out stage the system breaks up into hadrons 
and their resonances. The table of baryon and meson 
resonances implemented in the model is quite extensive. 
This allows for better accounting of the influence of final- 
state interactions on the generated spectra. The hard 
part of the model accounts for jet quenching effect, i.e., 
radiation and coUisional losses of partons traversing hot 
and dense media. The contribution of soft and hard pro- 
cesses to the total multiplicity of secondaries depends on 
both centrality of the collision and its energy and is tuned 
by model parameters to RHIC and LHC data. 

The paper is organized as follows. A brief description 
of the HYDJET++ is given in Sec. II. Section III presents 
the results of calculations of both V2 and V4 for charged 
particles in both considered reactions. The even compo- 
nents of the anisotropic flow and their ratio R = V4/V2 
are studied in the interval 10% < a jog^o < 50% in four 
centrality bins. In Sec. IV the interplay between jets and 
decays of resonances, as well as the roles of resonance de- 
cays in better realization and the jets in violation of the 
number-of-constituent-quark scaling are discussed. Con- 
clusions are drawn in Sec. V. 



II. THE HYDJET-h-l- EVENT GENERATOR 

The Monte Carlo event generator HYDJET-I-+ [16] was 
developed for fast but realistic simulation of hadron spec- 
tra in both central and non-central heavy-ion collisions 
at ultrarelativistic energies. It consists of two parts. The 
FASTMC [17, 18] event generator deals with the hydro- 
dynamic evolution of the fireball. Therefore, it describes 
the soft parts of particle spectra with the transverse mo- 
menta -pT < 2GcV/c. The hard processes are simulated 
by the HYDJET model [19] that propagates jets through 
hot and dense partonic medium. Both parts of the HY- 
DJET -1-4- generate particles independently. 

To allow for really fast generation of the spec- 
tra the FASTMC employs a parametrized hydrodynam- 
ics with Bjorken-like or Hubble-like freeze-out surface 
parametrization. Since at ultrarelativistic energies the 
particle densities at the stage of chemical freeze-out are 
quite high, a separation of the chemical and thermal 
freeze-out is also implemented. The mean number of par- 
ticipating nuclcons Npart at a given impact parameter h is 
calculated from the Glauber model of independent inelas- 
tic nucleon-nucleon collisions. After that the value of ef- 
fective volume of the fireball Vef / , that is directly propor- 
tional to Npart, is generated. When the effective volume 



of the source is known, the mean multiplicity of secon- 
daries produced at the spacelike freeze-out hypersurface 
is calculated. Parametrizations of the odd harmonics of 
the anisotropic flow are not implemented in the present 
version of HYDJET-I— |-, whereas the elliptic flow is gen- 
erated by means of the hydro-inspired parametrization 
that depends on momentum and spatial anisotropy of 
the emitting source. The model utilizes a very extensive 
table of ca. 360 baryon and meson resonances and their 
antistates together with the decay modes and branching 
ratios taken from the SHARE particle decay table [20]. 
After the proper tuning of the free parameters, the HYD- 
JET-|--|- simultaneously reproduces the main characteris- 
tics of heavy-ion collisions at RHIC and at LHC, such as 
hadron spectra and ratios, radial and elliptic flow, and 
femtoscopic momentum correlations. 

The multiple scattering of hard partons in the quark- 

gluon plasma (QGP) is generated by means of the HYD- 
JET model. This approach takes into account accumulat- 
ing energy loss, the gluon radiation, and coUisional loss, 
experienced by a parton traversing the QGP. The shad- 
owing effect [21] is implemented in the model as well. The 
PYQUEN routine [22] gcincratcs a single hard NN colli- 
sion. The simulation procedure includes the generation of 
the initial parton spectra with PYTHIA [23] and produc- 
tion vertexcs at a given impact parameter, rescattering- 
by-rescattering simulation of the parton path length in 
a dense medium, radiative and coUisional energy losses, 
and final hadronization for hard partons and in-medium 
emitted gluons according to the Lund string model [24]. 
Then, the full hard part of the event includes PYQUEN 
multi-jets generated around its mean value according to 
the binomial distribution. The mean number of jets pro- 
duced in A + A events is a product of the numbcir of 
binary NN sub-collisions at a given impact parameter 
and the integral cross section of the hard process in NN 
collisions with the minimal transverse momentum trans- 
fer. Further details of the model can be found in 
Refs. [16-19]. 

It is worth mentioning recent important modification 
of the HYDJET-I— t-. After the measurement of particle 
spectra in pp collisions at LHC it became clear that the 
set of model parameters employed by the PYTHIA 6.4 
version had to be tuned. Several modifications have 
been proposed [25, 26]. The application of standard 
PYTHIA 6.4 in the HYDJET-|-|- led to too early suppres- 
sion of elliptic flow of charged particles at intermediate 
transverse momenta in lead-lead collisions and, therefore, 
to the prediction of a weaker V2 [27, 28] compared to the 
data. Recently, the HYDJET-^-l- was modified [29] to 
implement the Pro-Q20 tunc of PYTHIA. In contrast to 
calculations of elliptic flow presented in [27, 28, 30], all 
simulations of Pb -|- Pb reactions at LHC energies in the 
present paper are performed with the upgraded HYD- 
JET-h-h. 
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FIG. 1; (Color online) Transverse momentum dependencies 
(triangles) of (a) V2 and (b) 114 of charged hadrons calcu- 
lated within the HYDJET++ for Au -|- Au collisions at 
= 200 j4 GeV at centrality a/a^co = 20 — 30%. Histograms 
show flow of directly produced particles in hydro-calculations 
(dashed lines), total hydrodynamic flow (solid lines), and flow 
produced by jets (dotted lines) . 



FIG. 3: (Color online) V2{pt) (full triangles) and V4{pt) (full 
circles) for charged particles in HYDJET-I— I- calculations of 
Au -I- Au collisions at ^/s — 200 A GeV at centrality a/agco (a) 
10 --20%, (b) 20 --30%, (c) 30 --40% and (d) 40- -50%, 
respectively. Dashed lines show hydrodynamic part of the 
calculations. Data from [10] are shown by open triangles {V2) 
and open squares {V4). 
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FIG. 2: (Color online) The same as Fig. 1 but for Pb -|- Pb 
coUisions at ^ = 2.76^ GeV. 



III. V2 AND V4, FROM HYDRODYNAMICS AND 
FROM JETS 

For the investigations of the second and the fourth 
flow harmonics, ca. 60 000 gold-gold and ca. 50 000 
lead-lead minimum bias collisions have been generated 
at Vs = 200^ GeV and = 2.76^ TeV, respectively. 
The transverse momentum dependencies of V2 and 174 ob- 
tained for the centralities 20—30% are shown in Fig. 1 for 
RHIC and in Fig. 2 for LHC energies. 

Together with the resulting distributions for V2{pt) 
and v^Ipt) we present separate contributions coming 
from (i) hadrons directly produced at the freeze-out hy- 
persurface in the hydrodynamic part, (ii) direct and sec- 



ondary hadrons created after the decays of resonances, 
and (iii) hadrons produced in the course of jet fragmen- 
tation. Recall briefly the main features of the V2{pt) 
behavior in IIYDJET-I--I-. The elliptic flow rises up to its 
maximum at intermediate px around 2.5— 3 GeV/ c and 
then rapidly drops. This falloff is observed in experi- 
mental data also. In the model its origin is traced to 
the interplay between the soft hydrolike processes and 
hard jets, as was studied in details in [27, 28]. The ideal 
hydrodynamics demonstrates continuous increase of the 
elliptic flow with rising transverse momentum. Because 
of the jet quenching the jets also develop an asimuthal 
anisotropy that increases with the pt too; however, this 
effect is quite weak and does not exceed few percent. The 
particle yield as a function of the transverse momentum 
drops more rapidly for hydroproduced hadrons than for 
hadrons from jets. Therefore, after a certain px thresh- 
old jet particles start to dominate the particle spectrum, 
thus leading to a weakening of the combined elliptic flow. 
A similar tendency is observed in Fig. 1 and Fig. 2 for 
the V4 also, but, because of the quite weak signal in the 
hydrodynamic part, the effect of the V4 falloff is not as 
pronounced as that of the elliptic flow. 

As shown in Fig. 1 decays of resonances can change 
the elliptic flow of directly produced hadrons with px < 
3 GeV/c by 1-2% at RHIC and by less than 1% at LHC; 
see Fig. 2. For the 174 the difference between the two his- 
tograms is negligible; i.e., resonance decays play a minor 
role for soft parts of both W2(pt) and W4(pt) distribu- 
tions. At Pt ~ 2.5 GeV/c jets come into play and change 
dramatically the shapes of the elliptic and hexadecapole 
flows. 

It is worth discussing here details concerning the de- 
termination of the flow components in the experiment 
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FIG. 4: (Color online) The same as Fig.3 but for Pb -f Pb 
collisions at -^s = 2.76A GeV. Experimental data are taken 
from [31]. 

and in the model. In the HYDJET++ simulations the 
elliptic flow is connected to the eccentricity of overlapped 
volume of colliding nuclei. No fluctuations in the loca- 
tion of nucleons within the overlapped zone are consid- 
ered. Therefore, the flow is determined with respect to 
the position of true reaction plane. The next even compo- 
nent, V4, is not parametrized in the present version of the 
model; i.e., the hexadecapole flow comes out here merely 
due to the elliptic flow. Thus, it should also be settled 
by the position of the true reaction plane. Because of 
the absence of the fluctuations and non-flow effects, the 
ratio V4/V2 obtained on an event-by-event basis equals 
that extracted by separate averaging of and f 2 over 
the whole simulated statistics. 

In the experiment the situation is more complex. For 
instance, in the standard event plane (EP) method the 
event flow vector Q„ for n-th harmonic is defined as (see 
[3] for details) 

Qn = iQn,x,Qn,v) = m COS (n(/).,), siu (n^j)^ 

= (Q„cos(n*„),Q„sin(n«'„))- (2) 

The quantities Wi and (f>i are the weight and the az- 
imuthal angle in the laboratory frame for the ith par- 
ticle, respectively. From Eq. (2) it follows that the event 
plane angle can be expressed via the arctan2 func- 
tion, which takes into account the signs of both vector 
components to place the angle in the correct quadrant, 

= arctan2((3„^j^, (5„^2;)/n . (3) 

The nth harmonic w„ of the anisotropic flow at given ra- 
pidity y, transverse momentum pT, and centrality a /ugeo 
is determined with respect to the angle 

v„{y,pT,cr/crgeo) = (cos [n(0j - '^n)]) (4) 



FIG. 5: (Color online) Ratio W4/(«2)^ vs. pr for charged 
particles in HYDJETH — h calculations of Au + Au collisions 

at ys = 200^ GeV at centrality a/o-geo (a) 10 20%, (b) 

20 - -30%, (c) 30 - -40% and (d) 40 - -50%, respectively. 
Full circles denote the hydro-|-jet calculations, open circles 
show only hydro-part, and open squares indicate the rescaled 
experimental data (see text for details). 

by averaging (. . .) over all particles in all measured 
events. It is easy to see that the event plane angle for the 
elliptic flow does not necessarily coincide with that for 
the hexadecapole flow ^1^4. To compare our model results 
with the experimental ones we need, therefore, the data 
where the fourth harmonic is extracted with respect to 
the \E'2 rather than the ^'4 event plane angle. 

To demonstrate the development of both V2 and V4 
at different centralities, we display the flow harmonics 
for charged particles in heavy-ion collisions at RHIC and 
LHC energies in Figs. 3 and 4, respectively. The exper- 
imental data by the PHENIX (RHIC) and the ALICE 
(LHC) Collaborations are plotted onto the simulations 
as well. One can see here that HYDJET-I--I- overesti- 
mates the elliptic flow of charged hadrons with trans- 
verse momenta 2GeV/c < pt < 4GeV/c in both reac- 
tions considered. This indicates that simplifled combina- 
tion of ideal hydrodynamics and jets is probably enough 
to simulate first two even harmonics of anisotropic flow 
at pt < 2 GeV/c, whereas at higher transverse momenta 
other mechanisms, e.g., coalescence, should be taken into 
account for better quantitative description of the flow be- 
havior. 

The elliptic flow produced by the jet hadrons with 
Pt ^ 2 GeV/c is almost zero. Because of the jet quench- 
ing, the flow increases to 3 5% with rising trans- 
verse momentum; however, the jets alone cannot provide 
strong flow signal, say V2 « 10%, even at LHC energies. 
Since the W4 created by jets is also very small, it would 
be instructive to study how the admixture of jet hadrons 
can alter the V4/V2 ratio. 

The ratio i? = i'4/i'| as a function of transverse mo- 
mentum is presented in Figs. 5 and 6 for four different 
centralities in Au -I- Au collisions at RHIC and in Pb 
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FIG. 6: (Color online) The same as Fig. 5 but for Pb + Pb 
collisions at = 2.76^ GeV. 



+ Pb collisions at LHC, respectively. The final result is 
compared here to the ratio obtained merely for hydro-like 
processes and to the experimental data. As was men- 
tioned in [11], the measured ratio should be noticeably 
larger than 0.5. There are event-by-event fluctuations 
that increase R even if both flow harmonics are deter- 
mined by means of the ^'2 event plane angle. The in- 
crease occurs because of the averaging of both V2 and 
V4 over the whole event sample before taking the ra- 
tio. These fluctuations are lacking in the HYDJET-I--I-; 
therefore, the data used for the comparison are prop- 
erly reduced. See [11, 12] for details. It is seen that 
parametrized hydrodynamics with the extended table of 
resonances already provides V4/V2 ~ 0.6, which is higher 
than the theoretical value of i? = 0.5. Jet particles in- 
crease this ratio further to value R w 0.65 at RHIC and 
i? ~ 0.7 at LHC. While the ratio R is insensitive to the 
transverse momentum at O.lGeV/c < pt < 3GeV/c, 
at higher pT it increases with rising transverse momen- 
tum both in model simulations and in the experiment, al- 
though the RHIC data favor a weaker dependence. Thor- 
ough study of this problem within the hydrodynamic 
model indicates [12] that neither the initial conditions 
nor the shear viscosity can be accounted for the rise of 
high-pT tail of the distribution. It looks like this rise can 
be attributed solely to jet phenomenon. 

At LHC energy the increase of R with rising transverse 
momentum aXpr >i GeV/c is quite distinct. The differ- 
ence between the model results and the data visible for 
semiperipheral collisions at 40% < a/ageo < 50% can be 
partly explained by the imperfect description of the el- 
liptic flow at pT > 2.5 GeV/c; see Fig. 4. Also, the STAR 

results concerning the V2 are about 15 20% higher than 

the PHENIX data, and the HYDJET-h-h model is tuned 
to averaged values provided by these two RHIC experi- 
ments. Nevertheless, the effect of hard processes is clear: 
The hydrodynamic part of the code yields rather flat ra- 
tio Vi/v2, whereas the jets provide the rise of the high-p^ 



FIG. 7: (Color online) Upper row: The KEt jriq dependence 
of elliptic flow for (a) direct hadrons, (b) hadrons produced 
both directly and from resonance decays, and (c) all hadrons 
produced in the HYDJET-I— I- model for An + An collisions 

at ^ = 200^ GeV with centrality 20 30%. Bottom row: 

The KET/uq dependence of the ratios (112 /n,) /(ff/3) for 
(d) direct hadrons, (e) direct hadrons plus hadrons from the 
decays, and (f) all hadrons. 

tail. 



IV. NUMBER-OF-CONSTITUENT-QUARK 
SCALING 

The number-of-constituent-quark (NCQ) scaling in the 
development of elliptic flow was first observed in Au -I- 
Au collisions at RHIC [13, 14]. If the elliptic flow, D2, 
and the transverse kinetic energy, KEt = ttit — rriQ, of 
any hadron species are divided by the number of con- 
stituent quarks, i.e., = 3 for a baryon and Uq — 2 
for a meson, then the scaling in V2{KEt) holds up until 
KEr/riq w 1 GeV [32]. The observation of the NCQ scal- 
ing seems to favor the idea of the elliptic flow formation 
already on a partonic level. For instance, as pointed out 
in [33] , the scaling is broken if hadrons are produced in 
the course of string fragmentation, whereas the process 
of quark coalescence leads to the scaling emergence. On 
the other hand, as was shown in Refs. [27, 28], the ful- 
fillment of the NCQ scaling at ultrarelativistic energies 
depends strongly on the interplay between the decays of 
resonances and jets. Note that the breaking of the NCQ 
scaling at LHC was observed experimentally in [15, 34]. 

To demonstrate the importance of both resonance de- 
cays and jets for the formation of NCQ scaling we plot the 
reduced functions V2 /nq{KET/nq) for several hadronic 
species obtained in HYDJET-(-+ simulations of heavy-ion 
cohisions at RHIC (Fig. 7) and at LHC (Fig. 8) energies 

in centrality bin 20 30%. These distributions are then 

also normalized to the flow of protons, vl^/uq : v\/?>, 
to see explicitly degree of the scaling fulfillment. The 
study is subdivided into three steps. The flow of hadrons 
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FIG. 8: (Color online) The same as Fig .7 but for Pb + Pb 
collisions at ^ = 2.76^ TeV. 

straight after the thermal freeze-out in hydrodynamic 
calculations is displayed in left windows. Central win- 
dows present this flow modified by the final state inter- 
actions, i.e., decays of resonances. Finally, right windows 
show the resulting flow of hadrons coming from all pro- 
cesses. 

At RHIC energy, it looks like at given centrality the di- 
rect pions, protons and kaons are produced already obey- 
ing the scaling within the 5 10% accuracy limit; see 

Figs. 7(a) and 7(d). The scaling holds also after decays 
of resonances as demonstrated in Figs. 7(b) and 7(e). 
Its fulfillment becomes slightly worse when hadrons from 
jets are taken into account; however, the NCQ scaling re- 
mains valid within 10% accuracy at least for three main 
hadron species. The situation is drastically changed for 
the collisions at LHC. Here spectra of directly produced 
particles do not possess any scaling properties, as one can 
see in Figs. 8(a) and 8(d). After final-state interactions 
the scaling conditions for hadrons in hydrodynamic sim- 
ulations are restored, as displayed in Figs. 8(b) and 8(e). 
Even mesons follow the unique trend. Why? Spectra 
of many light hadrons, especially pions and protons, are 
getting feed-down from heavy resonances, whereas the 
spectrum of </> remains unchanged. The resonance boost 
makes elliptic fiows of light hadrons harder. As a result, 
the NCQ scaling is fulfilled in a broad range of KET/uq 
in the hydro sector of the model. In contrast, hard pro- 



cesses cause significant distortions of particle spectra and 
lead to violation of the scaling conditions; see Figs. 8(c) 
and 8(f), in accordance with experimental observations 
[15, 34]. 



V. CONCLUSIONS 

Formation of elliptic V2 and hexadecapole W4 flows of 
hadrons in Au + An collisions at a/s = 200^ GeV and 
in Pb + Pb coUisions at = 2.7QA TeV is studied 
within the HYDJET-|--|- model. This model combines the 
parametrized hydrodynamics with hard processes (jets). 
Therefore, the main aim was to investigate the role of in- 
terplay between soft and hard processes for the develop- 
ment of flow. Several features have been observed. First, 
the jets are found to increase the ratio R = V4/V2 for 
both considered heavy-ion reactions. Second, jets lead 
to rise of the high-pr tail of the ratio R. Such a be- 
havior is observed experimentally but cannot be repro- 
duced by conventional hydro models relying on ideal or 
viscous hydrodynamics. Third, the resonance feed-down 
significantly enhances the fiow of light hadrons and mod- 
ifies their spectra toward the fulfillment of number-of- 
constituent-quark scaling. The flow of particles produced 
in jet fragmentation is quite weak, thus jets are working 
against the scaling. Due to interplay of resonance and 
jet contribution, the NCQ scaling works well only at cer- 
tain energies, where jets are not abundant. Because jet 
influence increases with rising collision energy, just ap- 
proximate NCQ scaling is observed at LHC despite the 
fact that the scaling holds for the pure hydrodynamic 
part of hadron spectra. At higher collision energies scal- 
ing performance should get worse. 
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